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ABSTRACT 


The complexity of a one-dimensional, fixed capacity delivery problem is examined. An 
instance of the problem is specified by n pairs of points (trips) on a straight line segment and a 
capacity k. nerd pair of points represents a source and destination of a package and the capacity 
is the limit on the allowable load of the delivery car. The car is required to travel back and forth 
on the line segment picking up packages at sources and delivering them to their corresponding 
destinations, in such a way that the car’s capacity never exceeds k. The delivery car may not 
unload a package temporarily. Therefore, once a source is visited the package remains part of the 
car’s cargo until its destination is reached. The focus here is on a further restriction, that each 
trip has the same direction in the line segment. It is required to find a route for the delivery car 


that completes all deliveries in a minimum distance. 


A polynomial time algorithm is presented for the capacity = 1 problem. The capacity > 3 


problem is shown to be NP-complete. 
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CHAPTER 1 


INTRODUCTION 


The boundary between those problems that can be solved in polynomial time and those that 
likely have no polynomial solution has been of interest to computing scientists since Cook’s paper 
[1972]. This thesis presents a problem called the delivery problem that has a polynomial solu- 
tion in its most restricted version. However, by relaxing only one parameter, it joins the class of 
problems for which the only known algorithms require a number of steps that grows exponentially 


with the size of the problem. It therefore becomes apparently intractable. 


1.1. Informal Description of the Delivery Problem 


The most general form of the deiivery problem can be thought of as one of optimizing the 
path of a delivery car which must haul a number of packages from given sources to given destina- 
tions without temporary unloading. The car has a finite capacity and each package uses up a 
specific amount of this capacity. As will be seen, this problem is too general to be interesting. 
Various restrictions can be placed on the probiem and will be discussed in Chapter 2. However 


the focus of this thesis is on the following very restricted version: 


n pairs of points on a straight line segment are given, each representing a source and desti- 
nation of a package. A delivery car is required to run back and forth on the line segment, 
picking up and delivering these packages with the restriction that it can carry no more than 
k ( ka fixed parameter greater than or equal to one) packages at one time. The car may not 
unload a package temporarily. Therefore, once a source is visited, the package remains part 
of the car’s cargo until its destination is reached. A further restriction is that all deliveries 
point in the same direction. That is, if (s,d) is a source-destination pair, then d>s. It is 


required to find a route for the car which completes all deliveries in the minimum distance. 
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1.1 Informal Description of the Delivery Problem 2 


Note that this is a significantly more restricted problem than the general one because: 
4: The car travels in only one dimension. 


ae All packages take up an equal capacity, 1, of the car’s capacity, k. 


3. All deliveries point in the same direction.! 


Note also that the car may start anywhere and end anywhere so long as it travels a minimum dis- 
tance. This restricted problem is called the one way, one-dimensional, k capacity delivery 


problem. It is abbreviated as /Way del-k. 


A summary of necessary definitions and concepts follows. A much more thorough treatment 
of the following material can be found in any source dealing with complexity of algorithms. This 


summary is loosely based on Garey and Johnson [1979]. 


1.2. Time Complexity of Algorithms 


The time complexity of an algorithm is a statement of its time requirements as a function of 
the input size. Suppose a problem, R, is solved by an algorithm, A, with (time) complexity f. 
This means that for every integer n, A solves every instance of R with input length n, in time no 
greater than f{n). Clearly the time required depends on both the measure and the computer 
model used. Similarly input length is dependent on the encoding scheme for the algorithm. The 
resulting time complexity functions, however, vary by no more than a polynomial function. 
These distinctions are immaterial in the theory of NP-completeness, which is concerned with iden- 


tifying broad categories of algorithm complexity. 


A function f:/-+I is Ofg) if and only if there exists a constant c such that | f(n)| < ¢|g(n)| 
for all n greater than some finite positive integer N. 

A polynomial time algorithm is one whose time complexity function fis O(p) for some poly- 
nomial function p. Thus there exists a constant c such that for every’ integer 
n, |f(n)| < c|p(n)| for all input of length n. 


This thesis will denote this direction as forward and represent it in figures as upward. 
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1.2 Time Complexity of Algorithms 2 


A problem for which no polynomial time algorithm exists is considered intractable. 


1.3. Decision and Optimization Problems 


Frequently two different versions of a problem can be considered, one requiring an optimiza- 
tion, the other requiring a decision. This is the case for the delivery problem of this thesis. The 


two versions are: 


Optimization problem 
Instance: n source-destination pairs and a capacity k. 


Problem: Find a feasible route with least length that completes aij deliveries. 


Decision problem 
Instance: n source-destination pairs, a capacity k, and a bound B. 
Question: Does there exist a feasible route of length less than or equal to B, that completes all 


deliveries? 


A solution is feasible if it meets the restrictions of the problem. In this case a route is feasible if 


the car never carries more than k packages at a time. 


The set of all instances of a decision problem, Ff, contains as a subset those instances which 


result in a yes decision. This subset is the set of yes-instances denoted by Ry. 


In order to be precise, the theory of NP-completeness is applied to decision problems only. 
However any algorithm that solves the optimization version of a problem, can be converted to one 
that solves the corresponding decision problem by just appending a comparison of the result with 
the bound B. Frequently the converse is also true. Using the decision algorithm, a binary search 
can be done on the possible outcomes to determine the optimal bound B’*. If the set of possible 
outcomes is small enough, this search takes no more than polynomial time. Next the optimal 
solution is constructed sequentially from partial solutions in polynomial time. This is done with a 


step by step extension of intermediate results until the final solution is found. 
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1.3 Decision and Optimization Problems 4 


1.4. Polynomial Reducibility 


A decision problem FR is reducible to a decision problem Q if there exists a constructive transfor- 
mation which converts any instance J of R into a corresponding instance J’ of Q such that Jis a 


yes-instance of F if and only if J’ is a yes-instance of Q. 


If the transformation can be constructed in time bounded by a polynomial function of the input 


size, then Ff is polynomially reducible to Q. 


1.5. NP-completeness and Intractability 


There exist many provably intractable problems. However so far, those that have been pro- 


ven intractable are either undecidable or intractable even using the non-deterministic model for 
computation.” There is another class of problems suspected of being intractable, but no proof yet 
exists. Finally, there are the provably tractable problems. These last two classes are dis- 
tinguished as follows. 


NP is the class of all decision problems which can be solved in polynomial time using a non- 


deterministic model of computation.? 


P is the class of all problems that can be solved in polynomial time using a deterministic model of 


computation. 


Since any algorithm for a deterministic Turing machine can be converted to one for a non- 
deterministic Turing machine by ignoring the extra power of the non-deterministic model, the 


class P is a subset of the class NP. 


A fundamental open problem of computer science today is whether or not P is a proper sub- 
set of NP. Although it is widely conjectured that P # NP, there is yet no problem known to be 


in NP and yet intractable. 


*Definitions and discussion of deterministic and non-deterministic models of computation can be found in Garey and 
Johnson [1979]. 
’Though these definitions and that of NP-completeness to follow are properly stated in the form of language recog- 


nition problems on deterministic and non-deterministic Turing machines, for the purposes of this thesis the equivalent 
characterization in terms of algorithms is sufficient. 
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1.5 NP-completeness and Intractability 5 


Instead, following the example of Cook, attention is focused on those problems that are 
among the hardest problems in NP. Cook proved [1971] that every problem in NP can be polyno- 
mially reduced to one particular problem in NP called Satisflability. Subsequently Karp [1972], 
and then many others, proved that a great number of decision problems shared this property with 


satisfiability. The importance of this class of problems has been recognized and given a name. 


A problem is NP-complete if it is in NP and every problem in NP can be polynomially reduced 


to it. 


Proving that a new problem is NP-complete is simplified by Cook’s pioneering work. Since 
polynomial reduction is transitive, it is sufficient to show two things. A problem, R, is NP- 


complete if: 
iS Risin NP. That is, there is a polynomial time non-deterministic algorithm that solves R. 


DS There is an NP-complete probiem that is polynomially reducible to R. 


1.6. Preview 


In Chapter 2 the general delivery problem is discussed and some variations are shown to be 
trivially NP-compiete. Other related problems are also cited. In Chapter 3 a polynomial algo- 
rithm or the one way, one dimensional, capacity one delivery problem, (1 Way de/-1), is 
developed and analyzed. For k>3, the one way, one dimensional, capacity k, delivery problem, 
(1 Way del-k), is NP-complete. Proof of this is the content of Chapter 4. Chapter 5 contains a 


summary of results, some conclusions and observations, and describes some related open problems. 
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CHAPTER 2 


RELATED PROBLEMS 


2.1. The General Delivery Problem 


The general delivery problem does not restrict the dimension of the space in which the 


trips lie. Nor is the weight of each package or the total car capacity restricted, though they are 


assumed to be positive integers.1 The definition of the decision problem version of the general 


delivery problem is: 


Decision general delivery 

Instance: n pairs of points (s;,d;) in Euclidean m-space, an n-tuple W=(w,, --- ,w,) associating a 
weight w; with each pair (s;,d;), a capacity k, and a bound B. 

Question: Is there a route for the car of length no more than B, such that it completes all 


deliveries and the total weight in the car at any one time is <k? 


The 1Way del-k problem is a subproblem of the general delivery problem. Other subprob- 


lems may be generated by changing the restrictions on the general problem. 


In light of the large degree of generality, it is not surprising that the general delivery prob- 
lem is NP-complete. In fact, the problem remains NP-complete when one or two of three different 


parameters are significantly restricted. Consider the following possible restrictions: 


‘Real numbers are not considered since it is intended to show that these problems are inherently NP-complete, not 
just difficult because of the complexity of computing with real numbers, No loss of generality results from eliminating 
rational numbers since a problem with rational weights and capacity can be converted to an equivalent integer problem 
by multiplying by the greatest common denominator. 
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2.1 The General Delivery Problem 7 


parameter unrestricted value restricted value 


dimension positive integer m 1 


weights w;=arbitrary integer << k w;=1 
for all z for all 7 
capacity positive integer & 1 


Except when the array of weights is arbitrary and the capacity restricted to 1, all combina- 
tions of these restrictions have meaning. However, the problem remains NP-complete unless ail 


three restrictions are imposed simultaneously. 


The notation delivery{mD,wW,kcap) will be used here to represent the delivery problem with 
parameters m (dimension), w (weight), and & (capacity) replaced from the table above. For exam- 
ple, delivery(1D,arbW ,kcap) is the one dimensional, arbitrary weights, k capacity delivery prob- 
lem. Using this notation, the 1Way del-k problem becomes delivery(1D,1W,kcap) with the addi- 


tional restriction that all deliveries are in the same direction. 


2.2. Some NP-complete Subproblems 


Some NP-complete results are immediate. A non-deterministic algorithm can solve the gen- 
eral delivery problem in polynomial time. After guessing a solution, it is only necessary to check 
that the proposed route is feasible, and to add up the length of the route to ensure it is less than 
or equal to B. Since both these last two operations can be done in linear time, the general 
delivery problem is in NP, and hence ali subproblems are as well. Therefore only the second 


requirement for NP-completeness, a constructive reduction, will be sketched. 


2.2.1. The Delivery(2D,1W,1cap) Problem 


Theorem 2.1 


Delivery(2D ,1W,1cap) is NP-complete. 


proof: 
The reduction is from the Euclidean Travelling Salesman (ETS) problem for the piane. This 


problem is NP-complete {[Papadimitriou,1977]. Let an instance, /, of ETS be specified by a 
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2.2.1 The Delivery(2D,1W,lcap) Problem 8 


collection, C, of n cities, a complete set, E, of intercity distances defined by the discrete Euclidean 
metric,” and a bound B. The transformation to an instance of delivery(2D,1W,lcap) consists of 
constructing a very short trip, (s;,d;), starting at each city c; in C. The idea is that, if the trips 
are short enough, they will have negligible effect on the length of a route either between the cities 


in ETS or through the deliveries in delivery(2D ,1W,1 cap). 


Construction of an instance, I’ of delivery(2D,1W,icap) corresponding to J in ETS is as fol- 
lows. At each city c; in ETS construct a source s, in delivery(2D,i1W.icap). In order to ensure 
integer values throughout, the edges in EF and the bound B are enlarged by multiplying the dis- 
tance function by 8n and the new trips are given unit length. Thus 4(s;,s;) = 3n8(c;,c;).° At 
each s,; construct a trip (s;,d;) of length one. Finally, set B’ = 3nB + 2n. This construction 
can be performed in polynomial time and has the property that J is a yes-instance of ETS if and 


only if J’ is a yes-instance of delivery(2D ,1W,\cap). 


For suppose there is a route (co, ae Ca) of length < B for J. Then the path 
(Sa, "++ 48) has length < 3nB. Hence, by the triangle inequality of the distance function, the 
route (Sada ae sea ta) has a length < 3nB + 2n for J'. Thus [ € ETSy implies 


I' € delivery(2D ,1W,1cap )y. Conversely, suppose J/' has a route of length 
< B' =3nB + 2n. Then because of the scale factor, / has a route of length 


) 
SnB+ on B+ sce But all edges have integer length. Therefore J has a route of length 


3n 3 


lA 


lA 


B. Thus I’ € delivery(2D,1W,icap)y implies J € ETSy.0 


2.2.2. The Delivery(1D,arbW,kcap) Problem 


Theorem 2.2 


Delivery(1D ,arbW,kcap ) is NP-complete. 


1 


2 A A 7 
? This metric is defined for the plane by (21-20) 24 (yi-Yo) 4 where [x] is defined as the smallest integer 
greater than or equal to Z. 


3 &z,y) is the distance from x to y. 
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2.2.2 The Delivery(1D,arbW,kcap) Problem 9 


proof: 
The reduction is from the NP-complete problem Partition [Karp,1972]. Let J be an instance of 
Partition specified by a set, A = {a),...,4,} of n positive integers.* The corresponding 


instance, I' of delivery(1D,arbW,kcap) contains n trips T;=(1,2) of unit length, each with the 


n 
same source and destination where the weight of trip 7; is 2a;. Capacity k 1s set to yy a; and the 
i=l 


bound to 3. Then there is a route of length 3 if and only if all deliveries can be made in two for- 


ward passes. This is possible if and only if there is a partition of A into 2 equal parts. 0 


The fact that delivery(1D ,1W,kcap) is also NP-complete is a consequence of a slightly more 


restricted problem, (the one way version), proved NP-complete in Chapter 4. 


2.3. The Polynomial Subproblem 


When the specified restrictions are placed on all three parameters, the 
delivery(1D,1W,1cap) problem results. A car travels back and forth on a straight line segment, 
making the trips from s,; to d; one at a time. The car is allowed to start anywhere and end any- 


where. The problem is to find a shortest route. 


Because the car’s capacity is one, after picking up a package at its source, it must travel 
directly to the corresponding destination. Therefore, for a specific instance, every route contains 
the same trips. Only the extra distance between trips can vary. Hence, minimizing this overhead 


is equivalent to minimizing the route length. 


Let A be an n by n array, where A, , is the distance from the i* destination to the 7" 


source, and A;;=oo. The problem, then, is the same as finding a single cycle permutation 


LO SS aa ie 


n 
(jist ipanonn, 1) ofk(lpidesn 1 )'suchitthat Ai, — max {A, 1) is minimal. ° 


4 The problem Partition asks if there is a subset A’ of A such that +s a; = uy Q;. 
a,cA a,c A-A’ 
n 
> Finding a single cycle permutation (§;=1,4,1,---,8,) of (1, ...,m) such that sy A; ; ii | is minimai yields a 
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minimum tour rather than a minimum route. 
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2.3. The Polynomial Subproblem 10 


This formulation highlights the similarity to the Travelling Salesman (route) problem. But 
contrary to the the claim in an unpublished manuscript by Chang [1976], the problems are not 
equivalent. The distance array for the Travelling Salesman problem may have arbitrary entries. 
In delivery(1D,1W,lcap), the matrix entries are more restricted than Chang assumed. The res- 
triction to one dimension implies that the rows of matrix A are not independent, Thus 


delivery(1D,1W,1cap) is a subproblem of the Travelling Salesman problem. 


Gilmore and Gomory [1964] solved a closely related problem. Their O(n log n) algorithm 
for sequencing a one state-variable machine, when applied to the input of delivery(1D ,1W,l1cap), 
produces the shortest tour for the delivery car, rather than the shortest route. If a fixed starting 
point, 6, and an endpoint, e, are specified, their algorithm could still be used by including the trip 
(e,5) in the set of trips, and finding the shortest tour. But there are only n*_n possible start and 
end point combinations in delivery(1D,1W,\cap) that could yield optimal routes. Therefore 
Gilmore’s algorithm could be applied n?-n times to yield the optimal route with arbitrary start 
and end points. Though this solution takes O(n"log n) and appears inefficient, it proves that 


delivery(1D,1W,1cap ) is indeed in P. 


An approach quite different from Gilmore and Gomory’s yields an O(n log n) solution pro- 
vided that, in addition to the current restrictions, all trips have the same direction. This algo- 
rithm is described in the next chapter. It is likely that a better algorithm than the O(n°log n) 


method mentioned above is possible for the two directional case as well. 


2.4. Other Related Problems 


Common features shared by various forms of the delivery problem and each of the Travel- 
ling Salesman, the Euclidean Travelling Salesman, and Partition problems have already been 


illustrated. Many other problems bear some degree of similarity to the delivery problem. 


An elevator problem discussed in Knuth’s Volume Three [1973], and solved by Karp, consid- 
ers the optimal transport of people between floors, using a single elevator. This problem resem- 


bles delivery(1D,1W,kcap) but has additional restrictions that admit a polynomial time solution 
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2.4 Other Related Problems 11 


unlike delivery(1D,1W,kcap). Though these extra properties result in quite a different problem, 
the technique used to calculate a lower bound resembles that used in Chapter 4 to facilitate proof 


of NP-completeness of the k capacity delivery problem. 


Wong, Liu, and Apter [1973] proved that a very simple algorithm transfers records on a 
rotating drum in optimal or nearly optimal time. Their problem is related to 1Way del-1. How- 
ever the drum scheduling problem operates on a circle instead of a line segment. In addition, the 
drum rotates constantly in one direction, whereas in 1 Way del-1 the delivery car can move in both 
directions. These crucial differences mean that their solution does not carry over to the delivery 


problem. 


A number of scheduling and sequencing problems have some superficial resemblance to vari- 
ous forms of the delivery problem. But those problems fail to capture the position information 
inherent in the sources and destinations of 1Way del-k. Also, they often include other restraints 
such as a partial ordering or multiple processors that are not analogous to parameters of the 


delivery problem. 
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CHAPTER 3 


THE ONE WAY, CAPACITY ONE DELIVERY PROBLEM 


In this chapter a polynomial solution for the one directional delivery problem where capa- 
city equals one is developed. The complexity of the algorithm is analyzed and is shown to be of 


optimal order. 


3.1. Initial Definitions and Terminology 
In order to make the problem and its solution precise the following definitions are used. 


A trip, T is an ordered pair (s,d) of integers. The functions sre(T) and dest(T) select the 
first and second coordinates of T respectively. Let lf = {7,,...,T7,} be a collection of trips. 
The function ls selects the trip T; € TF such that src(T,;) < sre(T,) for every T; ET. Thus, 
ls([) selects the trip in T with the least source. Similarly, gd(I’) selects the trip with the greatest 
destination. That is, gd([) = 7, if and only if dest(T,) > dest(T;) for every i. The expression 
top(T) is an abbreviation for dest(gd(I')). Similarly, bottom(T) is an abbreviation for src{ls(T)). 
The interval spanned by the trips in I’, boftom(I’) top(L), is denoted A(T). A set of tripsT is gap- 


less if, for any p € J([) there exists an 7 such that src(T;) < p < dest(T;). Let T contain n 


trips and a@ be an arrangement of the set {1,...,n}. A route #, is an ordering of the n trips 
in I, according towa. Thus Ka — (To, ile Bie Ts): The length, Len, of a route FR, is: 
n n—l 
Len({R,) = )) 14a, S01 + ))\4e,-%o,,,|- The overhead, Ov{R,) of a route, R, is: 
i=1 i=1 
n-1 


Ov(Ry) = di l4a,- %a,,,]- The pairs (dy.,8q,) are the links of the route Ry. Links are 


backward links if dy > Sa, and forward links otherwise. A subroute of I is an arrange- 


ment of a subset of the trips inT. 
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3.1 Initial Definitions and Terminology 13 


3.2. Problem Definition 


1 Way del-1 is defined by: 
Instance: a set I of n trips in a line segment PQ such that s; < d; for every T; ET. 
Problem: Find the arrangement y such that 

Len(R.) = minimum{Len(R,)|a@ is an arrangement of l,...,n}. 


For the remainder of this chapter, all trips are in one direction (forward). 


3.3. Results 


n n 
For all arrangements a, the first term of Len(R,), is ¥} | da. = (yy | = 95 (4 - s;). This is 
i=1 i=1 
just the sum of the lengths of the trips in I’ and is constant for ali routes. Therefore, to minimize 


the route length it is sufficient to find the route with minimum overhead. This is equivalent to 


minimizing the sum of the lengths of the links. 


3.3.1. Lower Bound on Overhead 


Let FR, be any route through T. Then: 


n-1 
Ov(R 4) <7 5 | do, oy | 


‘—1 


n-1 n-l 
D7 Fa, i Depa 


t—1 I 


n n 
= psi he yis t+ Sa, 


t=1 tai 


3) 4; - S)s; - top(T) + bottom(L) 


tot fd 


IV 


IV 


Denote this lower bound by LB(I'). Denote an optimal route through T by R*. Then LB(T) is a 


lower bound for Ov(R*) and any route R that has Ov(R) = LB(T) must be an optimal route. 


1Way del-1 requires that the minimum length be found over the set of all possible start and 
end points. However a similar lower bound is useful for a route with fixed start and end points. 


Let R', be aroute with fixed start and end points 6 and e respectively. Then: 
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I 


n-l 
| 5 oe | a ede. =e ath (tere 
i=1 


n n 
= Yid;- Mos; + b= €. 


jal i 


Ov(R! 4) 


Denote this lower bound by FLB(T). 


Notice that LB(I) is achieved by a route R, if and only if T,, = Is(), T, = gd(I), and 


Tn 
712 Sy, for every 7. That is, a route achieves the lower bound LB(I), if and only if it begins 


at the least source, ends at the greatest destination, and all links are backward links. Similarly, a 


route achieves lower bound FLB(L) if and only if it travels backward from 6 to s,,, from d, to e 


, 


and between trips. 


3.3.2. Staircases 

A staircase for [ is a subroute ST = (St,,...,St,,) of trips in IT such that 
St, = Is(T), St, = gd([), and sre(St;, :) < dest(St;) < dest(St;,,) for every « < m-1. 
Lemma L3.1 


If lis gapless, then a staircase for T can be constructed in time Of jT | }. 


proof: 
Let [ = {T,,...,T7,} be ordered so that src(T;) < src(T,,,). A staircase for this gapless set 


can be constructed as follows. 


CSS. Hee Jee? 
eae 2e: 
while St; A gd(I) 
do 
while dest(T;) < dest(St,) 
do 
esa ae 
od 
Je ce sts 
Siete FT; 
t+ 3+1 


od. 


Notice that in the inner loop there must exist an 2 such that dest(T;) > dest(St,) since 


top([) > dest(St;). Also the first trip, 7, such that desi(T) > dest(St;) must have 
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srce(T) < dest(St,;) since otherwise [ would not be gapless. Since 7 is incremented at every pass 
through the inner loop, the algorithm examines each trip in I once and therefore requires time 
OULD Teo 
Lemma L3.2 
Let I’ be a set of gapless trips together with a fixed start 5 and end e such that 5, e € I(I). 
If b > e, then there is a route R’ through TI such that Ov(R') achieves the lower bound 


FLB(I), and thus is optimal for fixed start and end. 


proof: 

Let ST = {St,,...,Stm} CT be a staircase for T. Let H be those trips T € ([ - ST) such 
that src(T) > 6. Similarly, let L be the set of trips T' € (T - ST) such that sre(T') < 6. 
Then [ = (SPF iL): Order H and L from greatest destination to least destination yielding 
Peer aan le (Lee aie Now seth, (flee aly Sic els Shndiiae Hye). 


Since 5 > e, R' has only backward links and therefore achieves lower bound FLB(T). 


This lemma yields an upper bound on the overhead of the optimal route R* through a gap- 
less set [. Let p be some point, p € /(['). By lemma L3.2 there is a route R, which starts and 


ends at p and achieves FLB(I). Therefore, 


Ov(R,) = FLB(L) 


Denote this upper bound for a gapless set T by UB(T). 


3.3.3. Additional Length of a Route 


For any route 2 through I’, Ov{R) can be calculated directly from the forward links in R 


and its start and end points. This is the content of the following lemma. 


Lemma L3.3 


Let R, be a route through I which starts at src(T, ) = s,, and ends at desi(T, ) = 4, . 
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Then 
Ov(R,) = LB(T) + 23)(sa,,, - da,) + (Sa, - dottom(L)) + (top(l) - a,) 
fl 


where flis the set of forward links. 


proof: 
Let b/ and flrepresent the set of backward links and forward links respectively. 


n-1 
Ov(R,) = DS | dy 3 ga 24 | 


'=1 


= Dilday tagyict 2 APe,,5 i4e,) 
bf fi 


= di(4a, a "Oe ) a dilda, . $a,,,) ne 207 (84... fe dy) 
ft fl 


bf 


n-l 
= dale, a9 $a;,,) at ila, i dy.) 


Didi Gey am ne + 807 tea laa, | dy) 


i—I t—1 fl 


= LB) + 2)i(8e,,,- da) + (8a, -— bottom(T)) + (top(P) - dy ).0 
fl 


Thus Ov(R,) is greater than LB(L) by twice the sum of its forward links plus the difference 


between its start and end points and the points bottom(I) and top (I) respectively. 


3.3.4. Improved Bounds 


If f is a set of trips with gaps, then the lower bound LB(I) can be improved. The following 
lemmas show that an optimal route through [ has forward links through the gaps in a set of trips. 
Lemma L3.4 

Suppose [ = PP -++T, where [; and [,;,, are sets of gapless trips arranged from 

lowest to highest sources and I(I';) Fj /(j.1) = 9. Let the gap between [; and I';,, be 


called G; and have size g,. That is, the distance from top(I ;} to bottom(I;,,) is gj. Then if 


yt 
R is a route through, Ov(R) > LB(T) + 23) yy. 


f=] 


proof: 
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Let the start of R be 6 and the end be e.’ For every h such that G, lies above 6, R must travel 
forward across G,, since it must complete all trips. Similarly, for every / such that G;, lies below 
e, R must travel forward across G;. Now consider the gaps that are both below 6 and above e. If 
there are no such gaps, then the result follows by lemma L3.3. Otherwise, 


(6 — bottom(T)) > YS} g; and (top([)- ¢) > S) g;. Therefore, using L3.3 
G,<b e<G, 


Ou(R) > LB(T) + 23) a + 2¥5 a + (6 - bottom(l)) + (top(T) - €) 
G>6b Ga=e 


t-1 
SA ies pes 
tI 


This improved lower bound for the overhead of a route through the trips in 
Reali ae with gaps Gai 4s, Gy 7h 1s: called IE BT). 

An upper bound for a route through a set of trips Il’, possibly with gaps, can now be formu- 
lated. Let’ =V\L)--- Ul: be as in lemma L3.4, and let ST; be a staircase for ;. Let R, be 
a route that starts at bottom(l), completes each staircase, S7T;, in order from bottom to top, and 
then completes the trips in; — ST; from top to bottom in order of decreasing sources finally end- 


ing in I’,. The only forward links for this route are across the gaps G;. So by lemma L3.3 
Ov(R) = LB) + 2)ifs0,,, - da,) + (top (FP) - bottom(T)) 
fl 


= LB(T) + 25%, + (top(I) — bottom{L)) 


n n t-1 
= yd; - SY) 5; + 2s ot 
f—! t=. i 


Therefore the optimal route must have overhead less than or equal to this length. This length is 


denoted by JUB(T). 


Lemma L3.5 
For any set of one-directional trips I’, there is a route, 2*, of optimal length which starts at 


a point lower than or equal to its ending point. 


"It is sufficient to assume 6 = ere(T,) and e = dest(T,) for some # and j, since otherwise R can be made shorter 
by eliminating this initial and final travel. 
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proof: 

Let R, be a route through [ with start 6 greater than end e. Assume I =Ty| J --- JI; as in 
lemma L3.4. 

Case 1: Suppose 5 and e are in the same interval /([,). Then by lemma _ L3.3, 


= dy) + (top([)-— e) + (6 - bottom(l)). But all gaps, G; above 


i+] 


Ov(R,) = LB(T) + 20(sq 
jl 


b or below e must have forward links through them. Therefore, in this case, all gaps have forward 
links. Thus 
Ov(R) > LB(T +25, gi + (top(T) — e) + (b — bottom(L)). 


Since b>e this is greater than JUB(I). Therefore R cannot be optimal. 


Case 2: Suppose bE/(F,) and e€/([;) where K>j7. By lemma L3.3 


Ov(Rq) = LB(T) + 2))(s0,,,- 4a,) + (top() - e) + (6 - bottom(l)). 


All gaps above 6, and all gaps below e must have forward links. Therefore 


Ov(R,) > LB(T) + 23) 9; + 2339; + (top([) - e) + (b — bottom(L)) 


i>k <y, 
= LB(T) + 23)9; + 2¥)9; + 2(6 - e) + (top([) — 5) + (e — bottom(L)) 
i>k i<j 
But 2(b - e) = => gi + 2 Ss I(U;) + 2(b — bottom(l,)) + 2(top(f,;) - e) 
=) tiar 


Therefore 


Ov(R,) > LB(T) + oF 4; + 2 5) I(P;) + 2(6 - bottom(T;)) + 2(top(T;) - e) 
=I J<ic<k 


+ (top([)- 6) + (e - bottom(L)). 
But a trip, R' can be constructed that begins at top([,) and ends at bottom(I,) and has over- 
head no greater than Ov(R,). First construct a staircase, ST; for each set [;. The route FR! 


proceeds as follows: 
i for i from j down to 1, do the trips in ([ — ST;) in order of decreasing sources. 


2: for i from 1 to j do the trips in ST; in order. 
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3. for 7 from 7+1 to k-1 do the trips in ; by completing ST; in order, and then completing 


([; — ST;) in order of decreasing sources. 
4, for 7 from k to ¢ do the trips in ST; in order. 
5. for from ¢ down to k do the trips in (['; — ST;) in order of decreasing sources. 


In steps 1 and 5, there are no forward links. In steps 2 and 4 the only forward links are 
across gaps. In step 3, for each 2, there is a forward link from the destination of the final trip in 
(C; - ST;) to top({;) which continues across gap G;. This forward length must be less than 
g + J(C;). Therefore, using lemma L3.3, the overhead of R’ can be bounded above. 


Ov(R') < LB(T) + 25g, + 2S) I(0,;) + (top (T) - bottom(l,)) 
iA Patek 


+ (top([,;) — bottom(T)). 
This is less than or equal to the overhead of route R,. Therefore R‘ is at least as short a route 


as JG. 


3.3.5. Additional Overhead of Subroutes 


Suppose S, = (So) Sacked ed, is a subroute of I containing exactly the trips in I’ CT. 


It is useful to determine how good the ordering of S, is assuming that these are the only trips to 
be covered. The additional overhead, AOv, of a subroute S, measures the amount of overhead 
accumulated by S, in excess of the lower bound for the set of trips in S,. It is defined by 
AOv(S,) = Ov(S,) - ILB(I'). Thus AOv(S,) = 0 if and only if S, starts at the /s(I’), ends at 
gd(I' ) and has backward links except across gaps. Otherwise 


AOv(Sq) = 2))(8e,,,- 4a,)-2 D) 91 + (sre(Sa,) — bottom(I’)) + (top(I') - dest(Sq_)). For a 
fl Gel’ 


complete route R,, AOv(R,) = 0 implies that R, achieves the lower bound and therefore must 


be optimal. 
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3.4. The Solution to 1Way del-1 


The optimal route for a set of trips [ = ri ihe UF: is found by determining a route, 
R* such that AOv(R“) is minimal. This proceeds in three steps. First, a subroute, S* is con- 
structed such that 1) AOv(S*) = 0, 2) S* contains /s([;) and gd(I;) for all 7, and 3) S* is maxi- 
mal in the sense that there is no subroute S’ containing a proper superset of the trips in S* that 
also has AOv(S') = 0. The second step is to insert the remaining trips into S* in such a way 
that minimum additional overhead is accumulated. The final step is to prove that the route thus 


constructed is optimal. 


3.4.1. First Part of the Solution 


Let PS Tile Leah: where each [; is a maximai gapless subset. Assume 


lT; | = m(i), that T is arranged from lowest to highest set, and that the trips in each subset 


(ie ee ...,7; ,.} are ordered by increasing sre(T; ). By lemma L3.1 there exists a staircase 


tm (1) 


ST; = (St, , Tueoe » Sti) for each set of trips [;. These staircases can be joined together from 


lowest to highest to yield a subroute S of T, 


Ss (Sty, t Fe » Sty Sta, hte Le epicucee ) I@in Sty, ..): The additional overhead of S is 


5(2) 
zero. 


t 
Let W be the set I’ — (_J ST; of trips not in any staircase, and let W,,..., Wy be the max- 


i=1 
imum gapless subsets of W arranged from lowest to highest. Some of these W, can be inserted 


into S yielding a longer subroute of I in such a way that the additional overhead remains zero. 


Lemma L3.6 
Let S = (S,,...,5,) be the subroute built from concatenating staircases together, and 
denote the interval, src(Sj+ 4) dest(S,), spanned by the link from S, to S;,, by I(link,). If 
there exists a j such that /(W;) () /(link;) A @ then W;, can be inserted into S without 


increasing the additional overhead. 
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proof: 


Suppose [(W;) C) I(link;) 49. Because of the construction of the staircases, the intersection 
must be with a backward link, and therefore dest(S,) > src(S;,,). But W; is gapless. So by 
lemma L3.2, there is a route R, through W, starting at b=sup(1(W,) I (ink,) ) and ending at 
e = inf(I(W,)A1(link,;)) with only backward links. Let S’ be the route that results from 


inserting R, between S; and S;,,. Since S" still starts at is([) and ends at gd(I') and has only 


backward links except across gaps, the additional overhead remains 0. 0 


Let S’ be the subroute built by inserting into S all such W; that have non-empty intersec- 
tion with some stairway link as in L3.6. If there exists W; up to W;,, all intersecting the same 
stairway link, link;, then each of them can be inserted between St, and St;., in order from 
highest to lowest set. Because of its construction, S*, clearly satisfies the first two of the three 
properties listed at the beginning of this section, that is: 1) AOv(S*)=0 and 2) S* contains 


Is([;) and gd(I’;) for every 2. 


3.4.2. Second Part of the Solution 


If every set W; has been inserted into S* then S* is a route for T with AOv({S*) = 0. It 


therefore has overhead that achieves JD B(I) and hence is an optimai route. 


In general however, there may remain some subsets which do not intersect a staircase link. 
Rename these uninserted sets from lowest to highest V,,...,V,, and denote the collection 
{V,,...,Vy} by A. For each of these sets the interval, /(V;,), is contained in the interval s; d;_ 
for exactly one trip S; in S*, and J(V,) intersects no other trips outside of V,. Figure 3.1 illus- 
trates the type of configuration possible for some of these remaining sets which are represented by 


ellipses. 


It is important to observe that the sets in A are determined by the structure of I’ and are 
independent of the staircase construction. A different staircase yields exactly the same set of trips 
as those in S*, and the same uninserted sets V,,...,V,. Also, as the following lemma shows, 


y) v 


S* satisfies the third property, namely 3) S* is a maximal subroute for which the properties 1} 
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Figure 3.1 Sample configuration of some V-sets 
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and 2) above hold. 


Lemma L3.7 
S* is a maximal subroute for [ such that 1) AOv(S*) =0 and 2) S* contains /s([;) and 


gd(T;) for every 7. 


proof: 

Let T be any trip in T - S* and let S; € S* be the trip such that /(T) C 1(S,). Let Vr be the 
V-set in A containing T. Consider any subroute, S' of I that includes both S$; and T and also 
includes /s([';) and gd(I';) for every 7. Then either S; precedes T or T precedes S,. 

Case I: Suppose S; precedes T in Ch bs 

Let T' be the first trip after T in S’ such that src(T'} > dest(T). If any such T' exists, then 
S' must contain a forward link at least as great as the distance from dest(T) to 
min {src(T'), dest(S;)}. If no such r! exists, then S’ must terminate at some point less than or 
equal to top(Vr). In either situation, AOv(S’) > 0. 

Case II: Suppose S; succeeds T in S’. 

Let T'' be the last trip before T in S’ such that dest(T’') < src(T). If any such T'' exists, 
then S’ must contain a forward link at least as great as the distance from 
max {dest(T'’), src(S;)} to src{T). If no such T’’ exists, then S' must begin at some point 
greater than or equal to bottom(Vr). Again, in either situation, AOv(S') > 0. 


Thus no trip can be added to those in S* while maintaining an additional overhead of 0. 0 


This part of the solution inserts the trips in A in such a way that they add a minimal 
amount of additional overhead. Each one of the sets in A is added to the subroute S* in one of 


three ways. 


i V, can be prepended to S*. That is no trips lower than bottom(V;) are completed before V; 


is begun. 


2 V, can be appended to S*. That is no trips higher than top(V,) are started after V, is 


completed. 
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3: V, can be inserted into another part of the subroute. That is either trips with sources lower 
than bottom(V,) occur both before and after V,, or trips with destinations higher than 


top(V,) occur both before and after V,. 


Methods 1. and 2. change the start and end points respectively but can be done without 
adding forward links. Method 3. does not change the start and end points of the subroute but 
adds forward links. Therefore any of the options add additional overhead to the new subroute. It 
is necessary to determine which of these three options is optimal for each V,, and construct a 
final route R* from this information. It must also be shown that any other route through T° 


would accumulate at least as much additional overhead as R *. 


The following lemma presents a way to prepend trips to S* {corresponding to method 1.) 


without need for any forward links, and calculates the additional overhead due to this construc- 


tion. 

Lemma L3.8 
The sets V,,...,V; can be prepended to S* yielding <i>S$* in such a way that 
AOv(<i>S*) = bottom(V;) - bottom(L). 

proof: 


Let ST(V;) be a staircase for V;. Let the subroute <i>" be defined by prepending the trips in 
sets V;,...,V; to S*. <i>S* begins at bottom(V;). It first completes the trips in ST(V;,), 
then those in V; - ST(V;) in order of decreasing sources. Next for & from i—1 to 1 it completes 
the set V; in order of decreasing sources. Finally S* is appended as is to this initial subroute. 


There are no new forward links in <i >S* so AOv(<i>S*) = bottom(V;) — bottom(l). 0 


The trips in sets V;,...,V, can be appended to S* in a way analogous to the prepending 
of lemma L3.8. Thus after completing S*, the sets V, down to V;,, are appended in order of 
decreasing sources. Finally the set V; is appended by completing V; — ST(V,) and then climbing 


the staircase ST(V,). The following lemma results. 
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Lemma L3.9 
The sets V,,...,V, can be appended to S* yielding S*</> in such a way that 


AOv(S* <I>) = top(l) — top(V,). 


Let <i>S*</> be the subroute with the trips in V,,...,V,; prepended to S* and 
(-1 
V,;,...,V, appended to S* as in lemmas L3.8 and L3.9. All trips except those in LJ Vy are 
ctl 


included in this subroute. These remaining trips can be inserted into <i>S*</> and the addi- 


tional overhead due to these insertions can be determined. 


It is convenient to rename each of the remaining sets V,;,71+1<4</-1 in a way that 


reflects which staircase trip, S$ 


;, contains J(V,). Suppose J(S;) contains 


TOV, ) ace 1(V,) 44+ lo¢<rsi-l butnot J(V,-;) nor I(V,.4) 


Then alternatively label 


Vas VP as) V9 filo y! ie where p(j) = r-g+1 indicating that [(V' #) C 3; 4; for 
1<t<p(J). 


Some further definitions are necessary. Refer to Figure 3.2 which illustrates some of the 
sets of trips embedded in one staircase trip. Let Oy, be the greatest destination less than 
bottom(V’,). Note that 6, is top(V’, ) if k>1 and 6, is the destination of some element in 
<i>S*<Il>. If there is no destination less than bottom(V',), (only the case when 
10, j= 1. F=1), then set O;. = bottom(l’). Similarly let Os, be the least source greater than 
top(V',). Then o, is bottom(V', ) if k<p(j) and 7;,,, 18 the source of some trip in 
<i>S*<i>. If there is no source greater than top(V’ a: (only possible when 
i=v+1, V',=V, and is contained in gd(I) ), then set o, = top{I). Let 5, be the distance 
from ftop(V',) to bottom(V', ) for 1<k<y(j)-1. Also let 6, be the distance from 


i} ; ; ‘ . ¢ t . . 
bottom(V’ , ) to 6, and let ae be the distance from top(V ie) toa; 


The following lemma inserts V;,,,...,V,. into <2: >S*</> using exactly u(7) forward 


links for each collection of sets V’ {par we y! Fie It uses the smallest available set of forward 
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Figure 3.2 Notation used for V-sets 
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links, and determines the additional overhead due to this construction. 


Lemma L3.10 


Suppose [(V' md sh SO V2 ma) are contained in S; and let 6;, 0, , and 6, be defined as 


above. Then the sets V’ fpocee VE j,,) can be inserted into <i>S*</> yielding a longer 


subroute <i>S/</> in such a way that 


(J) 
AOv( <i >Sj<l>) = A0v(<i>S° <i>) + 25/6, -2 max {6; }. 
1 =o O<k<p(y) 


proof: 


Let m(j) be the index such that eee = by G)" and let ST, be a staircase for V’, for 


1<k<p(j). The extended subroute <i>Sj</> has the trips in V’, up to V’ inserted 


Jm(y) 


into <i>S*<Il> after @;,, and has V'. up to V', inserted into <i >S‘</> before 


Im{s)+1 aly) 


Form a subroute S,, that first travels up the staircases ST; yore and then com- 


Feu) Jm(s) 


pletes the trips in V’, 


ST 56) down to V' ja ST;, in order of decreasing sources. Also con- 


mu) 


struct the subroute S;>5 that first completes the trips in Oy Ted ath down to 


(a laps pl; ae in order of decreasing sources, and then travels up the staircases from 
m(y)+ m(y 


ST, to ST; uy Then the subroute <:>S/</> is the subroute <i>S*</> with S;, 


Im(j)+1 


inserted after the trip in <i? >S*</> terminating at O;, and S,» inserted before the trip starting 


at a; (iy Notice that within each set V’ he there are no forward links. The only forward links are 
aU 


through the gaps oF for 0<k<m(j)-1, and m(j)+ 1<k<y(j). Thus 


u(y) 
AOv(<i>Sj<l>) = AOvu(<i>S*<i>)+ oS 6, -~2 max 6,;.0 
i 0<k<n(y) * 


Denote the route with all the sets V;,,,...,V;1 imserted into <i>S*</> using the 


method of lemma L3.10 as <i>S*"* <I>. 


Notice that for each j, the sets V;,..., Lae are inserted with y(j) forward links. The 


forward link due to each V, is either the distance (boitom(V,)-6,)= 6, or 
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(7; — top(V,)) = 6,. Let IC(V,) be the cost of inserting V,; as in L3.10. Then for any V, 
such that k<m/(j), where m(j) is the index of the maximum 6,, V, is inserted after 6, at cost 


IC(V;) = 26 For V, such that k > m(j), V; is inserted before o, at cost [C(V,) = 26,. 


Jey" 


Using this notation, 


A0v(<i>S*<l>) = (bottom(V;) — bottom(l)) + (top(U) — top(V,)) + y" IC(V;,). 


=t+1 
For every combination of 7 and / such that O0<?i<i/<v+1, there is a route <i>S* </> 
constructed using the methods of L3.8, L3.9, and L3.10. Let R’ be the one with least additional 


overhead. Thus 


AQuUhs)\— Janin p {ACvl<t > S <1 )} 


0<i<i<v+1 


= min | > IC(V,) + (bottom(V;)-bottom(l)) + (top (T)-top(V,))\ 


Pea a Perea i) 
R* can be constructed directly without building all routes of the form <i>S**</>. To 


determine those sets that accumulate less additional overhead by prepending to S* as opposed to 


j 
inserting, it is sufficient to find i such that (bottom(V;) - bottom(f)) - S$) IC(V,) is minimal. 
kl 


The procedure used to determine the optimal 2 is included in the next section where the compiete 


algorithm for finding R* is outlined. The idea is that the first index p such that 


p 
bottom(V,) — bottom(T) < })IC(V;) indicates that the sets V,,...,V, are less expensive to 
k=1 


prepend than to insert. After this point, further savings are detected by searching for a new point 


p’ 
p' > p such that (bottom(V,) — bottom{V,)}< 3) IC(V;,). Since the search is always for 
k=p+1 


further savings, the largest p’ found indicates those sets V,,..., V, that should not be inserted. 


Similarly, to determine those sets that accumulate less additional overhead through append- 


ing to S* than through inserting, it is sufficient to determine / such _ that 


Uv 
(top(') — top(V,)) - S)IC(V;,) is minimal. 
k=! 
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If /<z then, for the collection V;,...,V,;, prepending or appending both surpass inserting. 
Then R*.= <i>S*<i+1> for some i, and 
AOv(R*) = (bottom(V;) — bottom(l)) + (top (IT) — top(V;+1)) 
= (top(I') — bottom(T)) — (top(V;,1) — bottom(V;)) 
But this is least when (top(V;. 1) — bottom(V;)) is greatest. Therefore R* can still be easily deter- 


mined by finding 7 such that (top(V;,,) — bottom(V;)) 1s maximal. 


3.4.3. Optimality of PR’ 
It remains to show that f° is an optimal route through I. 


Theorem T3.1 


AOv(R*) = min{AOv(R)| Ris any route through I}. 


proof: 


Let Ry =(T,,.-.-,To,) be any route through I. Let 6 = srce(T,,), and e = dest(T, ). If 


1 


e <b then AOv(R,) > AOv(<i>S"*<i+1>) > AOv(R*). Therefore assume e > 6. It is 


sufficient to show that AOv(R,) > AOv(R*). The points b, and e partition the set Vy,...,V, 
into three parts: 1) V,,...,V,, where bottom(V,) < 6, 2) V,,...,V, where top(V,) > e, and 
3) Vii ---,V,-4. Consider the sets in the third partition. If there are none then 
AOv(R,) > (6 — bottom(l)) + (top(P) - e) 
> (bottom(V,) — bottom(L)) + (top(T} — top(V,)) 
= A0u(<q>S#*<r>) 
> AOv(R*) 
Otherwise R, requires at least one forward link for each of the sets V' ae Lig thptercis oe Vn ens 


since each J(V' a C S, for one and only one 3. But the shortest possible forward link for each 


V', is either (bottom(V' ,)-6,) or (7, — top(V’,). The route <qg>S*'<r> has the pro- 


perty that precisely the smallest possible forward links of these occur and no others. Therefore 


AOv(R,) > AOv(<q>S"%<r>) > AO(R*).0 
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3.5. Statement of the Algorithm 


The optimal route R* is constructed by building up a linked list of the trips in the correct 
order. The basic structure of the algorithm is outlined but the detail of lower level routines is 
omitted. The bookkeeping necessary to execute these procedures without backtracking is not 


included. 


Algorithm for 1Way del-1 


Step 1: 
order the trips yielding [ = (T,,...,T,) 
such that src(T;) < sre(T;, 4). 
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Step 2: 

# build S*, the maximal subroute such that AOv(S*) = 0 # 
subR + T, # put T, in the subroute # 

Ex + @ # Ex is the set of current extra trips # 
For 7 from 2 ton 

do 


if src(T;)<dest(subR ) 
then 
# in same gapless set # 
if dest(T;)>dest(subR ) 
then 
# a higher staircase trip # 
add_to _stair( 7) 


else 
# not a staircase trip # 
it Ez 
then 
# T,; becomes the current extra trip # 
start_new_Ex(T,) 
elif top(Ez) > src(T;) 
then 
# T; goes into current extra set # 
add_to_Ex(T;) 
else 
# T, does not intersect Ex. 
therefore new Ex set # 
# deal with current Ex set # 
if [(Ez)nlink(7)40 
then 
# Ex can be put into subR # 
insert(Ex,}) 
else 
# add Ex to list of V sets # 
add_to_V(Ex) 
fi 
start_new_Ex(T,} 
fi 
fi 
else 


# new gapless set # 
# deal with Ex set as above and # 


add_to _stair(T,) 
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3.5 Statement of the Algorithm 


Step 3: 

# Find those sets V,,...,V; that are better appended than inserted. 
This step is specified in more detail than others, 

as promised in the previous section. # 


# initialization # 
last_prepend_pt — bottom(I) # bottom of last set prepended # 


highest_prepend_set + 0 # index of last set prepended # 
total_IC + 0 # insert cost so far # 
1+] # index for current V-set # 
while < v # still V-sets left # 
do 
currST + stair(V;) # stair containing V;# 
maxgap «+ lowergap(V;) # gap below V,=6,_.# 
last_low_insert <— 72-1 #index of last set inserted from below# 


while I(V;) C currST 
do 
# decide whether to prepend or insert # 
# first find best insertion -- i.e. below or above # 


best_insert min {maxgap, uppergap(V,)} 
if (bottom(V;)-last_prepend_pt) < (total IC + 2 best_insert) 
then 
# prepend V; # 
highest_prepend_set + 7 
maxgap + uppergap( V;) 
last_prepend_pt +- bottom(V;) 
total_IC — 0 
else 
# better so far to insert V; # 
if maxgap > uppergap( V;) 
then 
# V;, best so far inserted before o( V;) # 
insert_pt (V;) —‘‘before of V;)’’ 
total IC + total IC + 2 uppergap( V;) 
else 
# V; and all lower sets in S; that are not 
prepended, are inserted after 6(V,) since new maxgap # 
for k from last_low_insert+ 1 to 2 
do 
# change insertion point for V.# 
insert_pt (V,) <— ‘‘after 6(V;)” 
od 
last_low_insert <— 7 
total_IC — total IC + 2 maxgap 
maxgap + uppergap (Vj) 
fi 
1 + 1+ 1 # ready to iook at next V; # 


od 
od 
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Step 4: 
# Find those sets that are better appended than inserted. 
This step has a structure mirroring Step 3. # 
Step 5: 
# If some sets are both below the highest prepend point and 
above the lowest append point then find the largest gap in 
the V-sets. # 
if highest_prepend_set > lowest_append_set 
then 
highest_prepend_set + 7 such that (top(V;, ,) - bottom(V,)) 
is maximum 
lowest_append_set + 74+ 1 
fi 


Step 6: 

# Build R* # 

prepend (Vi, vey Vhighest preena) 

append OV inaest append» +: °-» V5 

for k from highest_prepend_set+ 1 to lowest_append_set-1 


do insert (V;) 
od 


3.6. Time Complexity of the Algorithm 
The steps outlined in the algorithm can be examined for timing requirements. 


Step 1: 


Since this is just a sort by sources, time required is O(n log n). 


Step ic: 
Each trips is tested once in the order returned from Step 1. With the appropriate bookkeeping 


each trip can be dealt with in constant time. Therefore this step requires O(n). 


Step 3: 

Step 2 can pass all the required values for each V-set. The inner ‘‘while’’ loop increments 7 on 
each iteration. Since 2 is never decremented, the outer loop ensures there are no more than v 
passes through the inner loop ali together. The only part of this loop that includes more iteration 
is the ‘‘for’” loop. But each set can be examined here at most once. Therefore this step requires 


only O(v) time. 
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Step 4: 


As for Step 3, O(v) is sufficient. 


Step 5: 

A single pass through each V-set record locates the largest (top({V;,,)— bottom(V;)), and there- 
fore O(v) time is used in this step. 

Step 6: 

Staircases need to be built and pointers rearranged as previously determined. If bookkeeping is 


thorough, this requires only O(n) time at most. 


Since v < n, steps 2 through 6 execute in O(n). Thus ordering the trips is the most expen- 


sive part of the algorithm, which operates in time O(n log n). 


The following lemma shows that O(n log n) is a lower bound on the time required to solve 


the 1Way del-1 problem, and thus the algorithm is of optimal order. 


Lemma L3.11 


O(n log n) is the optimal order for solving 1 Way del-1. 


proof: 
Let A be any algorithm that solves {Way dei-1. Let S = {s;,...,8,} be a set of n integers. A 


can be used to sort the elements of S as follows. Construct a trip 7; = (s;,s;+ = for each ele- 


ment s; € S. Let J be an instance of {Way del-1 consisting of the collection of these trips. Since 


the trips in J do not overlap, the algorithm A will produce the route R = (T, 


ei 


1; ) such 
that src(T; ) < sre(T; ..) for every j. This yields an ordering of S such that 
Tiss iter cupid ie Thus A can be used to sort S. But sorting cannot be done in less than 


O(nlogn) time. Hence algorithm A takes at least O(n logn) time. The given algorithm for 


1 Way del-1 achieves this lower bound and therefore is of optimal order. Oo 
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CHAPTER 4 


THE ONE DIRECTIONAL, CAPACITY K > 2, DELIVERY PROBLEM 


Although the one directional delivery problem with capacity one can be solved in polyno- 
mial time, the problem is NP-complete for capacity greater than or equal to three. This is the 
content of this chapter. The problem for capacity two is so far unsolved and is discussed further 


in Chapter 5. 


4.1. Definitions 


A set of trips, I, has the same meaning here as in Chapter 3. Other definitions, with the 
exception of route, such as src, dest, ls, gd, top, and bottom are also used here as in Chapter 3. 
For capacity k>1 it is no longer true that a route is made up of trips joined together by links 
between them. Let I’ be a set of n one-directional trips, and let a@ be an arrangement of the set 
{1,...,2n}.. For arbitrary fixed capacity, &, a route, &,, is an ordering, R =r, ,..., ae 
of the 2n points {s,,...,5,,4,,...,4,} of sources and destinations of trips in T according to 


the arrangement a. In addition, it is necessary that 
1 For every 2, s; precedes d; and 


oe For every 7, the number of sources up to and including position 7 minus the number of desti- 


nations up to and including position 7 is less than or equai to k. 


These two conditions ensure that packages are picked up before they are delivered, and that the 
car is never carrying more than k packages at one time. Notice that this definition, for the case 


k = 1, is equivalent to the definition of a route in Chapter 3. 
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4.2. A Lower Bound for Route Length 


Because trips can be combined together when capacity is greater than one, the overhead as 
defined in Chapter 3 no longer has meaning. Still a lower bound is useful; but one is derived here 


for the total route length rather than for the overhead. The length, Len, of a route is 
2n-1 
Len(R,) = D5 [ames tare 
i=1 
The sources and destinations of trips in I’ partition the line segment from bottom(T) to 
top(I’) into no more than 2n-1 pieces. Let py,po,...,Pq, (¢<2n) be the points corresponding to 


sources and destinations in I’, taken in order from least to greatest. Let p(z) be the number of 


trips containing the line segment p; p;,,. Since & packages at most can be carried at one time, 


and since all trips are forward, there must be at least es — 1 passes through this line seg- 


ment to complete these deliveries. Therefore a distance of at least (p;,1— p;): b | - 7 must 


be travelled for each line segment. Summing these over successive segments yields the following 


lower bound. 


Lemma L4.1 
A lower bound on the length of a route R for the /Way del-k problem on a set of trips 


T-={7j,..., 2n} partitioning the line segment with points p;,...,p, is 


Ay ae Bes Di): oy 


k — (pq - i). 


tt 
If there is a route that starts at the least source, ends at the greatest destination, always 


travels forward with a full car, and backward with an empty car, then that route passes through 


Di Pi41 exactly 2 At —1 times, and achieves this lower bound. Also, since every interval 


P; Pi, Must be travelled at least 2 — 1 times, if a route travels any interval more times 


k 


than that, it cannot achieve the lower bound. 
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4.3. Problem Definition 


For capacity = 1, the optimization version of the one directional delivery problem was of 
interest. To prove NP-completeness of the capacity > 3 problem, the corresponding decision ver- 


sion is necessary. The capacity = k, one directional decision delivery problem is defined as: 


Decision 1Way del-k 
Instance: A set T of n trips, {7,,...,7,}, such that srce(7T;) < dest(T;) and a bound B > 0. 


Question: Does there exist a route, A, through I that has length less than or equal to B? 


4.4. 1Way del-k is NP-complete for k >3 


The reduction is from the Exact Cover by Three Sets, (XC), problem. This problem 


was proved NP-complete by Karp[1972]. XC is defined by: 


X3C 

Instance: A finite set X with |X| = 3g and a collection I] = {C,,...,C,} of three element 
subsets of X. 

Question: Does II contain an exact cover for X? That is, does there exist a subcollection II’ C I 


such that for every zt € X, z € C, E Il’ for exactly one ? 


Notice that IT’ is an exact cover for X if and only if LJ Ci=X and |IT' |=q. 
cel 


4.4.1. Proof of Completeness for Capacity Three 


Theorem T4.1 


1Way del-8 is NP-complete. 


proof: 
1Way del-k, and hence 1Way del-8, are in NP since the general delivery problem is in NP and this 


is a subproblem. 


To show that 1Way del-3 is complete, X3C is transformed to it. Let an arbitrary instance 


of X8C be specified by X = {z,,...,23,} and IIT={C,,...,Cn} such that each 
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C= {21 ,T1,,%,} where m € X. To construct an instance of 1Way dei—3 from this, a line seg- 


ment m + 3q + 1 units long is required, and points labelled at unit distances. One endpoint 
and the following m-1 points are consecutively labelled C,,C>o,...,C,. The following 3q points 
are labelled 2,22, ...,23,. The second last point is labelled out and the last point in. See Fig- 
ure 4.1. This line segment is assumed directed from C, to in, and trips are embedded in this line 
segment always in this direction. The set I’ of constructed trips consists of 6m members, which 


can be partitioned into 3 sets. 
The first set, 4, only depends on the size of X. For every i from 1 to 3q, construct the 
rips, — (zane eal == {7 eal ag 


The second set, [*, depends on the number of occurrences of each z, in the collection I. 
Let n(z;) be the number of times that z; appears in the collection II. For each 7 from 1 to 3g, 


construct 7(z;)-1 identical trips, 7,7; = (2;,out); 1 <j < n(z,)-1. Let? = eS ie ak 


three trips T;?; = (C;,2,) for j = 1,2,3. Let Fr? = U {T?,}. Now set T = Ws Uy te | Ba and 


1<li<m 
1<j<3 
let: 
m 
B = 2)S}(out - C,) + 2q(in — out) -(in - C,) 
[=1 
m 
= 2m-3q + 2951+ 2g -(m+ 3¢ + I) 
I=1 
=6mq + m?-q-1. 
The complete construction is illustrated in Figure 4.1. 
Note that |I''| = 34, and |f%| = 3m, and that each trip in these sets can be constructed 


in constant time. Since there are m 3-element subsets containing elements of X, the sum of the 


3g 
number of occurrences of alt z; in | JC, is 9m. Therefore, })(n(z;)-1) = 3(m-q) = |I*|. These 


1 


trips can also each be constructed in constant time. Determining 7{z,;) for all 7 requires a single 
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Figure 4.1 1lWay del-3 reduces to X3C 
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pass over the sets in II, and_ therefore takes time _ proportional to m. So 
|| = 3q¢ + 3m + 3(m-q) = 6m, and the construction of this instance of 1Way del—3 from an 
instance of X8C is polynomial in m. It remains to show that II contains an exact cover of X if 


and only if I has a route of length less than or equal to B. 


By lemma L4.1, a lower bound on the length of a route for [ is 


m ¢ 3q 
=2))1- Bk bi pote pelt ee -~(m + 3q + 1) 
— 3 i 3 3 


= ofeucn st, 3q:'2m + 2g - (m+ 3¢ + 1) 


=6mq + m*-q-1 
But this is the bound B. Therefore it is sufficient to show that [I contains an exact cover if and 
only if I has a route of length exactly B. 

Suppose that II’ CII is an exact cover of X. Let a route R for TF proceed as follows. R 
starts at C, and picks up all 3 trips T?,. If C, = {21,,71,,71,} € I! then as the car passes each 
destination of ‘eee it completes that trip and immediately picks up the trip (2.1) in I!. These 


trips are completed by travelling to in. Otherwise, if C, ¢ II’, then, as the car passes each desti- 


nation of T},, it immediately picks up one of the trips (2, ,out) in T°, These trips are all com- 


pleted by travelling to out. In either case, the car has started at C, and travelled directly to in or 
out always maintaining a full load without backtracking. If the car ends at zn after this first for- 


ward pass, it travels empty from in back to out. 


For each remaining C, that is not in II’, the car now travels from out to C; and back to 
out. From out to C, the car is empty. At C; it picks up all 3 trips T);, and changes direction 


heading for out. As it passes each of its current destinations, it picks up the trip (2; ,out) eT? 


and proceeds to out. 
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Now for each remaining C, € II’, the car moves directly back to Cy, and then forward to in. 
Travelling backward to C,, the car is empty. At each C, it picks up the three trips T; ; and 
changes direction. As it passes each of its current destinations, it immediately picks up the 


corresponding trips (7% ,2n) € I and proceeds to in. 


Since this route starts at the least source and ends at the greatest destination, travels for- 
ward with a full car, and backward with an empty car, it must achieve the lower bound. Thus 
the route has length B.' Therefore, if II contains an exact cover then there is a route for I with 


length less than or equal to B. 


Conversely, suppose there is a route, A, for [ with length B. Since B is the lower bound, R 


travels through each interval p; p;,,, exactly 2 a — 1 times, and therefore forward through 


each interval at most ali) times. But by the construction of I’, each interval p; p;,, has a 


multiple of 3 trips passing through it. Each interval C, C,,, is contained by 8S! trips; each inter- 
val from C,, to out is contained by $m trips; and the interval out in is contained by 9g@ trips. 
Therefore, when F travels forward through an interval, it must have a full car. Also, for a trip 7 


and any interval p; p;., C T, the package for T must be carried through that interval only one 
time. Otherwise all p(i) trips in p; pj, cannot be completed in Ae forward passes. This 


implies that when R travels backward, it must have an empty car. 


Now consider the $q trips containing out in. R must complete them in gq forward passes, 
each time executing 3 deliveries. On any one pass through out in, let T,,T,,Ty be the trips 
being executed. These trips have unique sources, say 2,,2,,2, since all trips to in have unique 
sources. But then the three trips completed by R just prior to T,,7,, and T, must have ter- 


minated at z,,Zy,, and z,, since otherwise & would have to travel forward without carrying full 


capacity or backward without being empty. Let these trips be T’,,7T',, and T',, and their 


1 This can also be seen by calculating the length of the described route directly. 
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sources be C,,C,, and Cy in {Cy,...,C,}. Since no trips terminate in this set, R must have 
been travelling backward before these trips. Hence C, = C, = Cy since otherwise again R 
would be travelling forward not full or backward not empty. Similarly the remaining g—1 passes 
through out in isolate a unique point C, € {Cy,..., Cn}. The sets in II corresponding to these 
points form a cover for X of size g, which must therefore be an exact cover. Thus if there is a 


route R for I of length less than or equal to B, then II contains an exact cover for X. 0 


4.4.2. Proof of Completeness for Capacity Greater Than 3 


It remains to show that {Way del-k for any k > 3 is also NP-complete. This can be done 


by supplementing the previous construction with additional trips. 


Corollary C4.1 


1Way del-k, for any k > 3, is NP-complete. 


proof: 
1Way del-k is clearly in NP since a non-deterministic algorithm need only guess a feasible 
arrangement of the 2n points consisting of sources and destinations of trips and check that the 


corresponding length is less than or equal to B. 


ESC is again used for the transformation. Let I be the set of trips constructed from an 
instance of ESC in the previous proof. For each /, 1</<m, define & —3 identical trips 
T',; =(C,out) for 1< 7 < k-3. Also define (k — 3)q identical trips T’’; = (out,in) for 
1<i < (k-3)q. Let If be the set of trips supplemented with these extra (k — 3)m + (k-3)q 
trips defined above. Set the bound B to be the same as defined in theorem T4.1. Then it is 
straightforward to see that there is an exact cover for X if and only if there is a /capacity route 
through I; of length B. The extra trips enforce a full car for all forward motion, and prevent 


simultaneous completion of trips starting at more than one point C;. 0 
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CHAPTER 5 


CONCLUSIONS AND FURTHER RESEARCH 


5.1. Observations Concerning the Delivery Problem 


The question of why a problem is NP-complete, could have direct bearing on whether 
P = NP. Therefore the answer will typically be difficult. However, there is one part of the 
1Way del—-k problem for a constant & that can be isolated as being polynomially solvable, thus 
indicating that this part alone is not responsible for the difficulty of the problem. Suppose some 
oracle declared the optimal order in which to pick up the sources of a given collection of trips. 
Because of the two restrictions on feasible routes, 1) sources precede corresponding destinations, 
and 2) no more than k deliveries can be current, there are well defined limits on the intermediate 
insertion of destinations into the route. A dynamic programming approach requiring only polyno- 
mial time can be used to place the destinations optimally. This is shown for capacity k=2 in the 
section of this chapter entitled ‘‘Open Problems’’. In this case the insertion of destinations has 
complexity O(n*). The generalization to arbitrary fixed capacity k results in an algorithm with 


ay aK Similarly, the sources can be inserted optimally into a fixed ordering of the 


complexity O(n 
destinations in polynomial time. It appears, therefore, that the problem’s difficulty is due to the 


way that sources and destinations interact and cannot be attributed to either part individually. 


Another observation involves the complicated solution to 1Way del—k presented in Chapter 
3. A lot of apparent complexity is due to seeking an optimal route rather than an optimal tour. 
If just a tour had been required, the solution would have been much simplier, because the extra 


details to accommodate special properties of starting points and ending points would disappear. 
uel n n t 

The solution becomes one of deriving the lower bound ))d; — })s; + 25) 9; for the overhead of 
i=% i=1 i=1 


any tour through I’. This is just a corollary of lemmas L3.2, and L3.4. Then it is easy to show 
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5.1 Observations Concerning the Delivery Problem 44 


that this lower bound is always attainable using the same method as used to derive the JUB in 
Chapter 3. This type of situation is not particularly rare. Frequently such a lack of symmetry 


hides a lot of inherent simplicity. 


5.2. Corollaries to the Completeness Result 


The NP-completeness of the two directional k capacity, delivery problem, 
delivery(1D,1W,kcap ) follows from the completeness of 1Way del—k, since the one way problem 
is a restricted version of the more general two way problem. Alternatively, this result can be pro- 
ven directly using a modification of the reduction constructed in Chapter 4. For capacity 3, it is 
sufficient to construct 3q return trips from in to out and 3 trips from out to each of the C,'s 
except C,. Then there is an exact cover if and only if there is a route such that the car is always 


full. The extension to capacity k simply mimics the extension for the one way case. 


Another variation, called the direct delivery problem, requires that once a package is 
picked up at s;, the car must travel directly to d; without any backtracking. It is allowed to pick 
up further packages before delivering 7; only if they are en route to d;. Even though this is a 
further restriction of 1Way del—k, the direct delivery problem remains NP-complete. In fact the 
proof in Chapter 4 that E3C reduces to 1Way del—3, actually constructs an instance of the direct 


delivery problem, thus proving the claim. 


5.3. Open Problems 


The problems for capacity equal to 2, lWay del-2 and the two way version, 
delivery(1D,1W,2cap), are unsolved. A completely different approach seems necessary in order to 
determine the status of these problems. The usefulness of E3C for attacking this problem 


appears exhausted at capacity equal to 3. 


However if the optimal ordering of the sources is known, then there is an algorithm that 
inserts the destinations optimally in O(n?) time. Let S = (s,,...,8,) be the oracle’s ordering 


of sources for the optimal route. Construct a multi-stage graph with n stages. Each node in the 
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5.3 Open Problems 45 
j" column (stage) represents a possible configuration in the car immediately after picking up the 
source s;. Thus there is one node for every possible set of undelivered destinations. At the Hy 
stage each of these sets contains destination d;. A node may also contain 1 other undelivered des- 


tinations selected from the j —1 trips already picked up. Therefore the 7” stage contains 


hoe (S)=i nodes. Figure 5.1 illustrates the multi-stage graph for n = 5 when capacity 


k = 2. Nodes may have a size of 1 {just d;) or 2 (a full car). Nodes of size 2 have three succes- 
sors since either 1 or 2 different destinations are visited before going to the next source. Nodes of 
size 1 have two successors since the current destination can be either visited or not before going 


to the next source. Therefore the number of edges going from nodes at stage 7 to nodes at stage 


j+lis2+4 3(7-1) = 37 - 1. 


{d,} —————+ {4} {d3} —______—_» { 44} 
{d,,do} {d,,d3} {dd} {d,,d5} 
{do,d3} {d>,d4} 


{d2,ds} 


{d3,d4} {d3,ds5} 


{d3,d5} 


Figure 5.1 Multi-stage graph for 5 trips 
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5.3 Open Problems 46 


Since an edge between a node at stage 7 and a node at stage 7+ 1 represents the delivery of 
0, 1, or 2 trips, the optimal ordering for these deliveries can be found in constant time, say K. 
Let the cost of calculating all edges between stage 2 and stage j be denoted C,; ;. Then the cost of 


calculating all edges between stage j and stage 7+ 1 is bounded by C, ;4,; < K-(3j — 1). Sum- 


n-1 
ming over all stages yields C,, < ))K-(3j - 1) = O(n’). A final stage of the multi-stage graph 
j=1 


is required to represent the delivery of the destinations remaining after s, is picked up. This cost 
is bounded by O(n). Thus the destinations can be inserted optimally into the ordering of the 


sources in time proportional to n’. 


Other variations could be considered. If the packages on the car waiting for delivery form a 
stack or a queue the problem changes significantly. It is not known if this restriction is sufficient 
to make the problem a member of P. As previously noted, the variation where the car travels on 
a simple closed curve, has been studied in the one directional case. It seems likely that the 
corresponding two directional case can be solved by combining some results from both the two 
directional line segment version and the one directional circular version. In any case, except when 
a particular application warrants it, pursuing numerous slight variations does not appear to be a 


useful endeavour. 


5.4. Applications 


The problem delivery(1D,1W,icap) originated as a manufacturing industry problem which 
was erroneously thought to be NP-complete. Efficient delivery systems are one application of the 
work presented in this thesis. However the primary results of this thesis, determining a set of res- 
trictions that can be applied to a problem without significantly simplifying it, are at present 


essentially of theoretical interest. 


Te 1 + ages 
af RED, oT 2) Yl gn oh ' 


or Fe "y> 
ean 


k Rate leah is ea ad 5 
_ 7 


pi at De rota Ly en 


VO) oe if = up A Ki = aa blory 
best, pe -s a 
+ eoosteniesh a9 to yrs vileh odd saoeenge 


ssotteniesh of! endT {nyo 


e - . 
“a af lena 


i etshmaor od biva 


‘ ciagis eogands areldeny Sd saree 


- aoe 
7 RJ 


to 


>} | a qt im ed ar maida a . 


= 


‘tule wed «ed orien beads, 


; _ 
5 -m va bh oe 30 Gn Fae | lancticonth wal gi 
_— 
a a 


mis?s7th 80 4 S hoe #o wIDY 1D ge sit bi 


— _— 


7 
A deonew aelonlgqe ial 


, £& t » foie to ( o> t, Kf Ji jyriaed fat Ay sf 


. s ——S~ - 
i Tek 1 n a4 el ol tava! ‘f mye Wed 
Lusnmnitg silt yaeewokt eind alata Gli 


d Jui age ige-Atre motcoig & o8 vest qa ott nad daa ' 


revistat Lesa 


7 


REFERENCES 


Chang, Robin [1976], “An overview of the Frito-Lay system,” unpublished manuscript. 


Cook, S. A. [1971], ‘The complexity of theorem proving procedures,’ Proceedings of the Third 
Annual ACM Symposium on Theory of Computing, Association for Computing Machinery, 
New York, 151-158. 


Garey, M. R., D. S. Johnson [1979], Computers and Intractability A Guide to the Theory of NP- 
Completeness, W. H. Freeman and Company, San Francisco. 


Gardner, Martin [1973], “Up and down elevator games and Piet Hein’s mechanical puzzles,” in 
Mathematical Games, Scientific American, February 1973,106-109. 


Gilmore, P. C., R. E. Gomory [1964], ‘Sequencing a one state-variable machine: a solvable case 
of the traveling salesman problem,” Operations Research 12, 655-679. 


Karp, R. M. [1972], ‘““Reducibility among combinatorial problems,”’ Complezity of Computer Com- 
putations, R. E. Miller and J. W. Thatcher (eds.), Plenum Press, New York, 85-103. 


Knuth, D. E. [1973], The Art of Computer Programming, Volume 8: Sorting and Searching, 
Addison-Wesley, Menlo Park, California. 


Liu, P. C. [1976], Analysis of Non-Planar Graphs, Ph.D. dissertation, Stevens Institute of Tech- 
nology, Castle Point, New Jersey. 


Papadimitriou, C. H. [1977], ‘‘The Euclidean traveling salesman problem is NP-complete,”’ 
Theoretical Computer Scrence 4, 237-244. 


Wong, C. K., C. L. Liu, J. Apter [1973], ‘‘A drum scheduling algorithm,” Proceedings of the l, 
Fachtagung uber Automatentheorie und Formale Sprachen, July, 1973, Bonn, 267-275. 


47 


iqnoeterm budetdedae “miajegt ved-orinT ol 1 : 


bn sgn eo) “ sotlearig pareeny instooils 


KL. odd of Chive) & wired: oe) ee etrsleqintd fetes} oontiioh 


sutp leniandsent pv arski aot % bar, . riturg etgvals vor) bos qu" J 


“iL rand TG! oiel pofiete?. ei 4h hae oMoodisiinmvatak sahi's 


wnt 
naanie | a . 


saiteqmoD +) eoveloneAé yetogee’y 1 


7) 


ovetveetT ae? "ergo Sae aamostd HW y 


OLY) OTC! vanaded waahonek 4 —— eons? fe 


-@ 


stiducm oldsiver-oiete a0 ¢,poheasepye” [FORT sane 


wat 66 OI dapsone eraliseie sualdeng nemuentee g 


We Us ierry 2. ,m hiaid ds JNK CE OS, Prints qld 
dioY wa 22919 ceed Lbs) inivtedT .W 1. bos 9 45 


gerinot .% smiwloY enimivesgesS welaqne uk Oke bas 
sf wall ? 44 a laalft valeoWs 


> 


jien! evravoled molinsqantly Chel mow! tans aa ¥ ri ste ‘ 
tat all 


ol 
perm! wo.4 


alder uni get nadie? oft Sq or Se 


pC-Ths .b ewric? weleqereD Jlesty wath 


hiaaat irdtitogis gcrtubod i arid ie Ter] sedqg fh i Pini as i a Be 


PON MG heh ste 
whet nie bie 
Em ais, 


Ve 
ty 


MONA enleg 
Ste BANAT 
ue 


‘si 


ay; 


hy 


hia 
a4 


Ont 4 i 
age LenS Meg 
NG Sake Ait 


rials hye 


E 


#9 


ie 
in 


baat 


i a 
Hew 


crith 


waa 


wet 
Has Ling 
i 7 
Ps 


yee 


i 


ca 


aa 
Ten ean tReet oe : : wee 
re nea gee REL ga Motes ? 


: ie 
uy oe 
FP dee nite 
stanky nish 
AMF adhere 


Ha 


fat ‘ 
ge egy) at ty i n pe ae 
$ aH Hs Ra ‘ Setdybaned wpa) ot 
Piha bab en a oi Hrselye ins 


4, 
reer Nate lp 
‘ ae 


the q bee 


a 


a 


i 


nee nyu 

PR I OR tea 

Hatch ramp ban nap ee Veh tang 
Reet aL etic eps 


eit 
Shift oem 
AG De 
a 

My 


iD 
atsavich 


We 


whontpian 
i 


aM at 
hich 
WNcbpare 


ee a hu MPN, 
HANES WI AC oo 
a ; aud fad 


4 


it 
Cp RNA 
a 


7 
Ns Mea NT 
ia 


OM 
poke 


He ect 
See 


reas ine 
ba 


Thats 
4 


boom Ledy bate eg 
ah ut trae) 
Webiets 


a 
vor 


pores tn 


